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Abstract 

We study a branching Brownian motion Z in among obstacles scattered 
according to a Poisson random measure with a radially decaying intensity. 
Obstacles are balls with constant radius and each one works as a trap for 
the whole motion when hit by a particle. Considering a general offspring 
distribution, we derive the decay rate of the annealed probability that none 
of the particles of Z hits a trap, asymptotically in time t. This proves to be 
a rich problem motivating the proof of a more general result about the speed 
of branching Brownian motion conditioned on non-extinction. We provide 
an appropriate skeleton decomposition for the underlying Galton-Watson 
process when supercritical and show through a non-trivial comparison that 
the doomed particles do not contribute to the asymptotic decay rate. 
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1. Introduction 


Branching Brownian motion (BBM) among random obstacles function¬ 
ing as traps has been studied recently in [3, E] • We study a BBM 

that evolves in where a radially decaying field of Poissonian traps are 
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present. We are mainly interested in the trap-avoiding probability of the 
system up to time t, where the underlying Galton-Watson process (GWP) 
has a general offspring distribution. Investigation of this problem leads us 
to finding the speed of BBM conditioned on non-extinction, which is of 
independent interest. 

Let Z = [Z{t))t>o be a d-dimensional BBM with branching rate /3 > 0 
and offspring distribution A. The process starts with a single particle at 
the origin, which performs a Brownian motion in for a random time 
which is distributed exponentially with constant parameter j3. Then, the 
particle dies and simultaneously gives birth to a random number of particles 
distributed according to the offspring distribution A, which is a probability 
measure on N = {0,1,2,...}. Similarly, each offspring particle repeats the 
same procedure independently of all others, starting from the position of 
her parent. In this way, one obtains a measure-valued Markov process Z = 
{Z{t))t>o, where Zt can be identihed as a particle conhguration for given 
t>0. By assumption, Z(0) = 5o- The total mass process \Z\ = {\Z{t)\)t>o is 
a continuous time branching process with rate /?, and the number of particles 
in generation n oi \Z\ is a GWP (A^(n))neN with offspring distribution A. 
The initial particle present at f = 0 constitutes the 0*^ generation, the 
offspring of the initial particle constitute the 1** generation, and so forth. 
We denote the extinction time of the process \Z\ by r, which is formally 
dehned as r = inf {f > 0 : |Z(f)| = 0}, where we use the convention that 
inf0 = oo. We then denote the event of extinction of the process \Z\ by T, 
and formally write = {r < ooj. We use the term non-extinction for the 
event Let P be the probability corresponding to the process Z, and E 
the corresponding expectation. 

The branching Brownian motion is assumed to live in a random environ¬ 
ment consisting of Poissonian traps. Let be the d-dimensional Borel 

sets, and let 11 denote the Poisson random measure on with a spa¬ 

tially dependent locally finite mean measure u such that dv/dx exists and 
is continuous on and 


dx 


I 


\d-i ■ 


\x\ 


oo, I > 0, 


( 1 ) 


where dx is for the Lebesgue measure. A random trap configuration K with 
radius a on is defined as 

K= IJ B{xi,a), 

3;iesupp(n) 


where B{xi,a) is the closed ball of radius a > 0 centered at x*. Let P be 
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the probability for the Poisson random measure, and E the corresponding 
expectation. 

For A G and t > 0, notice that |supp(Z(t)) n A\ is the number of 

particles located in A at time t. For t >0, let 

= U supp(Z(s)) (2) 

se[o,t] 

be the range of Z up to time t. Now let T be the first time that Z hits a 
trap, 

T = inf {t > 0 : |supp(Z(t)) n iCl > 0} = inf {t > 0 : R{t) n iF / 0} . 


Then, the event of trap-avoiding up to time t ol Z among the Poissonian 
traps is given by {T > t}. 

We aim to analyze (ExP)(r > t), the annealed (averaged) trap-avoiding 
probability of the ^stem up to time t, for a BBM with a general offspring 
distribution A. In [7|, the asymptotic decay of the annealed probability has 
been found as t —>• oo in the strictly dyadic branching case, that is, A(2) = 1. 
On the way to its generalization, we find out that the problem proves to be 
rich enough to require several stand-alone results that we also contribute in 
this paper. The most important of these is on the speed of BBM conditioned 
on non-extinction as given in Theorem [TJ When the GWP is supercritical, 
the particles are grouped into those with infinite or finite line of descent, 
so-called “skeleton” and “doomed” particles, respectively. In other words, 
a skeleton decomposition is performed to analyze the problem. We show 
by a comparison argument that the doomed particles do not contribute to 
the conditional speed. Here, the speed of BBM refers to the growth rate of 
the radius of the minimal ball containing the range of the BBM. It has first 
been studied by McKean 14, 1^, later by Bramson 0,1 and Chauvin and 
Rouault 0], and is still subject to active research. 

In Theorem [21 we prove the asymptotic decay of the annealed trap¬ 
avoiding probability as a large deviation result. It is an important and non¬ 
trivial application of Theorem [H The problem of trap-avoiding asymptotics 
for BBM among Poissonian traps has been previously studied by Englander 
0 for the case d > 2 and where the trap intensity was uniform. Then, 
in search for extending the result to the case d = 1, Englander and den 
Hollander 0 considered the more interesting case where the trap intensity 
was radially decaying as given in ([1]) . Both of these works were for a strictly 
dyadic underlying GWP, and later Oz and Qaglar [l6| studied the uniform 
Poisson intensity case for GWPs with a general offspring distribution. In 
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the present work, we consider a radially decaying trap intensity as in (fTj). 
and a general offspring distribution for the underlying GWP. We show that, 
conditioned on non-extinction, the doomed particles of the decomposition 
of the supercritical GWP do not contribute at all to the asymptotic decay 
rate of the trap-avoiding probability. In comparison with the strictly dyadic 
case, the mean and another functional of the offspring distribution appear 
as parameters in the rate expression (see ([5])). We refer the reader to [8| for 
an interes ting survey article on the topic of BBM among Poissonian traps, 
and to [3, 12] for various related problems. 

Theorem [3] is a technical result identifying the rate function of Theorem [2] 
through tedious analysis. 

The organization of the paper is as follows. In Section 2, we give the 
statements of the three theorems described above as our main results. Sec¬ 
tion 3 includes the essential lemmas together with their proofs. In sections 
4, 5 and 6, we prove theorems 1, 2 and 3, respectively. 


2. Main Results 

To formulate our main results, we introduce further notation. Let / be 
the probability generating function (p.g.f.) of the offspring distribution, and 
be the mean number of offspring. 


CX) 


j=0 


and define 


j=0 

m = /.i — 1. 


Throughout this work, we assume that fi < oo and without loss of generality 
that A(l) = 0 (see the proof of Lemma H]). Now let q = P{£) be the 
probability of extinction for the underlying GWP, and let 


a = l- f\q). 


Note that A(0) =0 implies that q = 0 since two or more offspring is produced 
each time a particle branches. When A(0) = 0, it then follow that a = 1 as 
/'(O) = A(l), and A(l) = 0 by assumption. Also, q = I when // < 1, and 
q e (0,1) when A(0) > 0 and ^ > 1. 
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For r,b > 0, define 


gdir,b) = j ilfid-i ’ 

JBrO) + ber ^ 

where e = (1, 0,... , 0) is the unit vector in the direction of the first coordi¬ 
nate. Finally, let 

= = (3) 

Sd V 2m 

where Sd is the surface area of the d-dimensional unit ball (si = 2, S 2 = 2tt, 
S 3 = 47r, etc.). 

Finally, recall from ([2]) that Usg[o,t]^(s) denotes the range of Z up to 
time t, and define 

M(t) = inf {r > 0 : R{t) C 5^(0)} for d > 1, (4) 

to be the radius of the minimal ball containing R{t). Now we state our main 
results. 

Theorem 1 (Conditional speed of BBM). Suppose that the underlying GWP 
is supercritical, i.e., m > 0. Then, ford > 1, conditioned on non-extinction, 
converges to y/2Prn in P-probability as t ^ oo. 

Regarding the speed of BBM, we refer the reader to 0,0 for the case 
d = 1, and where A(0) = 0. We note that almost sure speed results exist 
too (see e.g. 0) , but for our purposes, convergence in probability suffices. 
For d > 2 and A(0) = 0, we derive the speed in Proposition [1] in Section 
4. Theorem [H which is also proved in Section 4, covers the general case 
where A(0) can be nonzero. In this result, we see that conditioned on non¬ 
extinction of the underlying GWP, the speed remains as y/2flm, which can 
be explained as follows. When the supercritical GWP is decomposed into 
skeleton and doomed particles (see Lemma [4]), the doomed particles do not 
contribute to the speed, and the skeleton particles alone form another GWP 
with the same mean as the original one. 

Theorem 2 (Variational formula). Fix d,f ,/3,a. Define 

min \ fiar] + ^ lgd{\/2fim{l - t/),c)| . (5) 

r,e[o,i],ce[o,y2;g] I 2r] J 


(For T] = 0 put c = 0 and gd{-\/2j3m, 0) = SdV2firn.) 
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1. If A(0) = 0, then 

lim ^log(E X P){T >t) = -I{l,f,l3,d). 

t—>-oo t 

2. If A(0) > 0 and m > 0, then 

lim i log(E X P) (T > t\n = -1(1, f, /3, d). 

t^OO t 

3. If A(0) > 0 and m < 0, then 

lim (E X P) (T > t) = (E X P) {T > t) > 0, 

t^OO 

where r = inf {f > 0 : |Z(f)| = 0} is the extinction time and £ is the 
event of extinction for the underlying GWP. 

In comparison with the variational result given in [3, Theorem 1.1] for 
dyadic branching, Theorem [2] contributes two extra factors, both naturally 
depending only on /, in the expression for the rate function ([5]). The extra 
factor a in the first term appears when A(0) > 0 (if A(0) = 0, then a = 1). 
It is a consequence of the modified p.g.f. for the offspring distribution of the 
skeleton particles, which are part of the decomposition of the supercritical 
GWP conditioned on non-extinction. The extra factor m inside gd appears, 
because the speed of BBM is now no longer -v/2^ but ^/2prn. Since a = 
1 = m for strictly dyadic branching, the result in is recovered by (l5|) as 
a special case. 

The next result is purely analytic, and solves the variational problem for 
the rate function /(Z,/,/3, d). We generalize the corresponding result in [^. 

Theorem 3 (Crossover). Fix f ,(3,a. 

1. For d > 1 and all I Icr, the variational problem has a unique pair of 
minimizers, denoted by rf = p*(I, f, (3, d) and c* = c*{l, f, f3,d). 


2. For d = 1, 


and 


l<hr: /(/,/,/3,d)=/3-^, 

^cr 

l > Icr ■ 1(1, f, (3, d) = 13a, 

I < Icr ■ t]* = 0, C* = 0, 
l>lcr: V* = 1, c* = 0, 


( 6 ) 


( 7 ) 
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where Icr = 



3. For d>2, 


and 


where Icr 


where 


l<lcr: I{lJ,l3,d)=/3^, 

^cr 

l>lcr- I{l,f,P,d) < 


I < Icr ■ V* = 0) C* = 0, 

I > Icr '■ 0 < rj* < 1, 0 < c* < 

Ad*cr! 


Id 


— 1 + ^1 + ImaM^ 
2mMj 


€ ( 0 , a), 


( 8 ) 


(9) 


( 10 ) 


= max [ 5 rf(i?, 0 ) - 1)]. (11) 

2Srf771_Rg(0,oo) 

Furthermore, c* > y/2l3m{l — rf) when I > Rr- 

Comparing this theorem to the corresponding one in [^, we see that in 
this case the extra factors a and m appear ubiquitously in the formulas, 
.^ain, we note that the formulas above reduce to the corresponding ones in 
[Tl] upon setting a = 1, m = 1. 


3. Preparations 

In this section we prove seven preparatory lemmas. Lemma [T] is needed 
directly in the proof of Theorem [2j Lemma [2] and Lemma [3] are needed in 
the proof of Proposition [Hand Theorem [H which are central in the proof of 
Theorem [2j Lemma 0] is needed in the proof of both Theorem [1] and Theo¬ 
rem [21 Lemma [5] and Lemma [6] are needed to prove Theorem [H Lemma [7| 
is used in the proof of Proposition [H 

Lemma 1. Let r > 0 and b > 0. Then, 
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Proof. By substituting y = x/t, we see that 


9 d{rt, bt) = 

Je 


dx 


Brt{0) \x + bte\ 


d-l 


= t 


dy 


'b.(o) \y + be[ 


JZT ='t9d{r,b). (12) 


Recall that integrals of asymptotically equivalent positive continuous func¬ 
tions are asymptotically equivalent if one of the integrals diverges. Apply 
this result to the radial integral. By assumption, dv/dx is continuous on 
As the d-dimensional volume element has radial component |x|'^“^d|x|, 
the integrands of both radial integrals below are continuous on M+. Hence, 
we have 

v{Brt{bte))= / / , dx = lgd{rt,bt). (13) 

JBrtibte) dx jBrdhte) 

Combining (fT^ and (fT^ yields lim ~ '''^hich implies the result. 

□ 

Lemma 2 (Overproduction). Let 5 > 0. Then for any t 

P{\Z{t)\ > < e-^\ 

Proof. Use the fact that E\Z{t)\ = (see [I^) and then apply Markov 
inequality. □ 

Lemma 3 (Comparison). Suppose that A(0) = 0. Let B C be open or 
closed. Let x ^ B. Let Px be the law of Brownian motion, denoted by W, 
starting from x, and let P^^ be the law of BBM starting from 6x. Define the 
first exit times from B 


V'B = inf {t > 0 : lU(t) G B^} 

ifB =inf {t > 0 : \supp{Z{t)) C > 1} . 


Then for any /c G N and t >0 

PsA^B > t\\Z{t)\ <k)> [Pxif^B > t)f. 

Proof. See the corresponding proof in , which is written for strictly dyadic 
branching, and note that the proof can easily be extended to the case of 
general A-GWP with A(0) = 0, since the population does not decrease for 
such processes, that is, if for any k G N+, \Z{t)\ < k, then \Z{s)\ < k for all 
s G [0, t]. CH 
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Lemma 4 (Decomposition of the supercritical GWP). Let N = (A^(n))neN 
be a supereritical GWP with offspring distribution X, where A(0) > 0. Let f 
be the corresponding p.g.f.. Let q be the probability of extinction for N, and 
define q = 1 — q. Let N*(n) be the number of particles in generation n that 
have an infinite line of descent. 

1. The law of N given extinction is the same as that of a GWP with p.g.f. 

f{s) := f{qs)/q. (14) 

2. The law of N* = (N* {n))n£'N given non-extinction is the same as that 
of a GWP with p.g.f. 


f{s) ■=[f{q + qs) -q]/q. (15) 


3. The law of N given non-extinction is the same as that of a process N 
generated as follows: Let N* be the process consisting of the particles 
of infinite line of descent. To each particle of N* having no children, 
add ni more children, that each start an independent GWP with p.g.f. 
f, where ni has the p.g.f. 


, , , _ {D^°f){qs) 
• (Dnoffiq) 


(16) 


where D is the differentiation operator, and all ni and all GWPs added 
are mutually independent given N*. The resultant process is N. 

4. Let X = (X(t))tgK+ be a continuous time branching process with rate 
fi such that the number of particles in generation n of X is N(n) 
for n G N. Let f be the p.g.f. for the offspring distribution of N, 
and let X* (t) be the number of particles with infinite line of descent 
present at time t. Then, the process X* = {X*{t))^^]^+ given non¬ 
extinction is equal in law to a continuous time branching process with 
rate j3{l — f'{q)) and offspring distribution, whose p.g.f. is 


f-.= [r{s)-r{q)s]/{i-nq)). (i?) 


Remark. The reader should note that in m, positive mass on 1 is allowed. 
For example, if X{0) = 1/3, A(2) = 2/3, f{s) = 1/3 + (2/3)s^, then q = 
1/2 = q, and f'{q) = 2/3, while f*{s) = (2/3)s + (l/3)s^. This distribution 
still has the same mean as originally (i.e., 4/3/, but here we put 2/3 weight 
on 1. 

What |77| ) says is that, equivalently (without putting mass on 1), it has 
rate /3/3 and p.g.f. f = s'^. 
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Proof. See [3, Proposition 4.10] for a proof of parts 1,2 and 3, and j^, 
Section 1.12] for further details on the decomposition of supercritical GWPs. 
For a proof of part 4, observe that a continuous time branching process with 
rate /3 and offspring distribution A is equal in law to a continuous time 
branching process with rate /3(1 — A(l)) and offspring distribution A, where 
A(l) = 0 and X{j) = XHVil — A(l)) for j / 1. To complete the proof, see 
the proof of part 3 of [iH, Theorem 1]. □ 

Lemma 5. If s € [0,1], then f*{s) < f{s). 

Proof. By (fT4]) and (fTSl) . we need to show that 

[f{q + qs) -q\/q< f{qs)/q. (18) 

The left-hand side of (fTS|l is equal to [f{s + q{l — s)) — q\/{l — q). Since / 
is convex, by Jensen’s inequality, it is enough to show that 

[s + (1 - s)f{q) - q]/{l -q)< f{qs)/q. 

Since / is greater than or equal to the identity on [0,g], s G [0,1] implies 
that f{qs) > qs. We now have, using f{q) = q, that 

s + {1 - s)f{q) - q ^ s+(l-s)g-g ^ ^ ^ f{qs) 


Lemma 6. Let Z\ = {Zi{t))t>Q and Z 2 = {Z 2 {t))t>o be two BBMs with 
the same constant branching rate f3 > 0, and p.g.f. ’s (pi and (p 2 for the 
offspring distributions, respectively. Let P^^ and be the corresponding 
probabilities for processes that start at time s with a single particle at position 
X G M'’*. Assume that (pi < (p 2 on [0,1]. Then, for any Borel set B, and any 
fixed time t > s, 

PUR\t) CB)< PlffR\t) c B), (19) 

where Rfft) is the range of process Zj up to time t for j = 1,2. 

Proof. Clearly, li x ^ B, then Pf ffR^it) C P) = 0 = Pf ffR^{t) C B), so 
the inequality in (I19p is trivially satisfied. Now, assume that x G P. Define 

u^{s,x,t) := PiffRfft) C P), s <t, for j = 1,2. 

Let N, j = 1,2, be the first branching time after s, for a single particle 
located at x G M'’* at time s. Then, is exponentially distributed with 
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parameter /3 for j = 1, 2. Since they induce the same corresponding 

sub-probability measure Q on [s,t]. Let E^^^ and E^^^ be the expectations 
corresponding to and Pg ^^ respectively. Then, conditioning on D and 
using the strong Markov property along with the branching property, we 
obtain for j = 1,2, that 

u^{s,x,t)= J 7r{r — s,x, B)K^^l(j)j (u^{r,Yr,t)) Q{dr) 

+ s,x,B), 


where 7r{t — s,x,B) is the probability that a single Brownian particle Y 
that starts at time s and position x G remains in B until time t, while 

(r) 

denotes the expectation conditioned on the event that Y remains in B 
within the time interval [s,r]. Thus, 4>i < 4>2 implies that 


{v} — u^){s,x,t) < / TT{r — s,x,B)x 

J S 

[01 {u^{r,YrR)) -01 {v?{r,YrR))] Q{dr). (20) 


Set s = 0 and use that Q has an explicitly known bounded density g = g{t) 
to arrive at 

-u‘^){0,x,t) < [ 7r{r)g{r)E‘^f^l [0i {u^{r,Yr,t)) - 0i {u^{r,Yr,t))] dr. 

Jo 

Letting n-^(x, t) := u^{0,x,t), j = 1,2, and making the substitution r i-)- t—r, 

pt 

{v^-v‘^){x,t) < / 7r(t-r)5r(t-r)Eo“'’^ [0i {v^{Yt-r,r)) - 0i {v^(Yt-r, r))] dr. 

Jo 

(21) 

Now let w := (n^ — n^) V 0, where V denotes maximum. Recall that 0i is 
monotone nondecreasing and convex on [0,1], having a bounded derivative. 
Hence 0i(n^) — 0i(n^) > 0 holds if and only if — v'^ > 0, in which case 
01 (n^) — 01 (n^) < C{v^ — n^) for some C > 0. Therefore, 

pt 

w{x,t)<C TT{t - r)g{t - r)E^Q~'^'^w{Yt-r,r)dr. (22) 

Jo 

Note that if w{x,t) = 0, then (I22p holds since the right-hand side is non¬ 
negative, and if w{x,t) > 0, then ([2^ holds by the dehnition of w and by 
(|^ . Define F{r) := EQ^^'"'^w{Yt-r,r), and note that F{t) = w{x,t). We 
have ^ 

F{t)<C / Tr{t — r)g{t — r)F{r) dr, 

Jo 
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and, by Gronwall’s inequality, we conclude that F(t) < 0. Hence, w{x, t) < 
0. But w{x,t) > 0 by definition, therefore w{x,t) = 0, that is — v'^ > 0, 
which means that < u'^, as claimed. □ 


Lemma 7 (Brownian hitting times in M'^). Let B be a Brownian motion in 
starting at the origin, with corresponding probability Pq, and let k > 0. 
Then, as t ^ oo. 


Po 


I sup |Hs| > kt 

\o<s<t 


exp{—^(l + o(l))} . 


(23) 


Furthermore, for d = 1, even the following stronger statement is true. Let 
mt = iniQ<s<t Bs o-nd Mt = supq< 5 <^Hs. (The process Mt — mt is the 
range process of B.) Then, as t ^ oo. 


Kt 

Po{Mt -mt> kt) = exp ——{1 + o(l)) 


(24) 


Proof. For d = 1, recall that [l0|, equation 7.3.3] for any a > 0, 


Po ( sup \Bs 

\0<s<t 



Setting a = kt above, we see via rHopital’s rule that 


^0 


I sup |Ps| > kt 

\o<s<t 


= exp<{ —^(l + o(l)) 


We now extend this result to d > 2. Note that the lower bound holds 
trivially since the projection of a d-dimensional Brownian motion onto the 
hrst coordinate axis is a one-dimensional Brownian motion. In order to 
prove the upper bound, we first prove ()24p . 

Let d = 1. Define 


According to p.l99 in 
6c satisfies 


6c = inf {s > 0\Ms — ms = c} . 

and references therein, the Laplace transform of 


E exp ( — 


A2 


2 "J 1-1- cosh(Ac) ’ 
Hence, by exponential Markov inequality, 


A > 0. 


A^ 


t 1 E exp ( ——6c 1 = exp 


Po{6c <t) < exp 


1 -|- cosh(Ac) 
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Taking c = kt, one obtains 


Po{Mt -mt> kt) = Po{9kt <t)< exp 



2 

1 + cosli{Xkt) 


Optimizing, we see that the estimate is the sharpest when A = A;, in which 
case, we arrive at the desired upper bound. Since Pq (supq< 5 <j ^ kt) < 
Po{AIt — mt> kt), the lower bound coincides with the upper bound, which 
implies the result. 

For d> 2, let us consider, more generally, the Brownian motion starting 
from X, denoted by B^, x £ and let r = |x|. We will use the well known 
fact [171, Example 8.4.1] that if R = \B\ and r > 0, then the process R, 
called the d-dimensional Bessel process is the strong solution of the one¬ 
dimensional stochastic differential equation on (0, oo): 

d - I 

Rt = r+ / ^—ds + Wt, 

Jo 

where W is the standard Brownian motion. 

Using the strong Markov property of Brownian motion, applied at the 
first hitting time of the p-sphere, it is clear that in order to verify the lemma, 
it suffices to prove it when B is replaced by B^, and |x| = p > 0. This 
follows, since for every fixed p such that 0 < p < kt, we have 


Po ( sup I Bs 1 > It ) < Pq ( sup I Bs I > kt] = Pp ( sup Rs > kt 

\0<S<t ) \0<S<t+T J \o<s<t 

(25) 

where r is the hrst hitting time of B on the p-sphere, and {Pr',r > 0} are 
the probabilities for R (or for W). 

Now consider B^ and define the sequence of stopping times 0 = tq < 
<7o < Ti < (Ti,... with respect to the filtration generated by R (and thus, 
also with respect to the one generated by W) as follows: 

To := 0; fJo := inf {s > 0\Rs = p/2} , 



and for i > 1, 


Tj+i := inf {s > a/Rs = p} ; fTj+i := inf {s > Ti+i\Rs = p/2} 
Note that for i > 0 and s £ [rj,(7j], 
d-1 


Bs — P + 


/ 

J Ti 


2R, 


dz + Ws-Wr,<p + p-\d - 1)AI + W, - Wr„ (26) 
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where A| := s — Ti. 

Since -Rs < /o for Uj < s < Tj+i and i > 0, it is also clear that for t > p/k, 
the relation supo<s<t Rs > kt is tantamount to 

sup sup Rs > kt. 

2>0 TiAt<s<aiAt 

Putting this together with (I25|) . it follows that Pp(supQ<^<j Rg > kt) can be 
upper estimated by 

Pp >0,3 Ti At < s < ai At : Wg — > kt — p~^{d — 1)A| — p) . 

This can be further upper estimated by 

Pp (^Mt -mt>[k- -j- - = ^0 -mt>[k- -J- “ > 

for any <5 > 0, as long as t > p/6, where the last equality follows from the 
translation invariance of Brownian motion. To complete the proof, hx p, 
<5 > 0, let t —)• oo, and use the already proven relation (1241) . along with (I25p . 
Finally let 5 ^ 0 and p —^ oo. □ 


4. Speed of BBM 

In this section, we first state a well-known convergence in probability 
result regarding the speed of BBM, which holds when the dimension is one, 
and when each particle is certain to give at least one offspring, i.e., A(0) = 0. 
We then extend this result to any dimension, and hnally to the case where 
A(0) > 0. The results of this section are central in proving Theorem [2j 
Recall from dH) that M{t) = inf {r > 0 : R{t) C Rr-(O)} for d > 1 denotes 
the radius of the minimal ball containing the range of the process up to time 
t, where R{t) stands for the range up to time t. 

Proposition 1 (Speed of BBM). Suppose that A(0) = 0. For d> I, 
converges to yJ2j3m in P-probability as t —>• oo. 


Remark. It is easy to see that the quantity fim is invariant under changing 
the setting from the ‘canonical one’ (that is, putting zero mass on 1) to a 
non-canonical one (that is, assigning positive mass to 1). 


Proof. For d = 1, the reader is referred to 14, [H 3)11(1 [^ • Now lot d 2. 
Since the projection of Z onto the 1st coordinate axis is a one-dimensional 
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BBM with branching rate /3, the lower estimate for convergence in proba¬ 
bility follows from the result for d = 1 and the inequality 



Ve > 0 and Vt, where P* denotes the law of the one-dimensional projection 
of Z. To prove the upper estimate for convergence in probability, let e > 0 
and let B = 5(y2^+e)f (0). Pick 5 > 0 such that 

+ e)^ > /3m + 6. (27) 

Recall that iI^b and tpB are the first exit times from B for a single Brownian 
particle starting at the origin and for a BBM, respectively. See that these 
events are identical: > y/2l3m + e| = < t|. Estimate 

+ < P{\Z{t)\ > 

+ P{^Pb < t\ \Z{t)\ < (28) 

By Lemma [21 the first term on the right-hand side tends to zero exponen¬ 
tially fast. Now consider the second term. Recall that Px is the law of 
Brownian motion starting from x, and Ps^ is the law of BBM starting from 
5x- By Lemma [HI we have the estimate 

P5o(^b > t\\Z{t)\ < [e(^-+^)*J) > [Po(V^B > 

By Lemma [71 


|Po(fe > = 


1 — exp — 


[(V2/3m -b £)t]" 
2t 


exp(o(t)) 


|^g(^m+5)tj 


(29) 

By (|27p . using the binomial expansion, right-hand side of (I29p tends to 1 
exponentially fast as t —)• oo, so that (|28P yields the desired upper estimate 


P + e 


0 , 


which completes the proof. 


□ 
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4 . 1 . Proof of Theoreml^ 

We emphasize that in this proof, the branching rate for aii types of 
particies for aii processes considered is constant and is equai to /3, and 
aii probabiiities written shouid be understood as conditioned on the non¬ 
extinction of the underiying GWP. 

First, decompose the supercriticai GWP into two sets of particies: par- 
ticies with infinite iine of descent, which, foiiowing the terminoiogy in [181 ] . 
we caii skeieton particles; and particles with finite line of descent, which we 
call the doomed particles. By Lemma U each skeleton particle gives rise to 
a GWP consisting only of skeleton particles with p.g.f. 

f*is) := [f{q +qs) - q]/q, (30) 

where q is the probability of extinction for the underlying GWP and q = 
1 — q. Also, by Lemma 01 each doomed particle starts her own single-type 
GWP with p.g.f. 

f{s) := f{qs)/q. (31) 

Observe that (/*)'(!) = f'{l) = /r. Let k = f'{l) = f'{q), and define 
k = K — 1. It is clear that each skeleton particle produces at least one skele¬ 
ton particle, and can produce doomed particles as well, whereas a doomed 
particle only produces doomed particles. 

Now we translate this decomposition into the language of multi-type 
GWPs, see [ 1 , Ghapter 5]. We define the skeleton particles to be of type 1, 
and the doomed particles to be of type 2, so that the underlying GWP is a 
two-type GWP, and we have the following decomposition for the BBM: 

{Zit))t>o = {Z\t),Z\t))t>o, 

where Z^ is the process consisting of the skeleton particles, and Z^ is the 
one consisting of the doomed particles. Note that Z^ by itself is a BBM, 
however, is not. Let \Z^(t)\, j = 1,2 be the number of particles of 
type j existing at time t. Observe that conditioning Z on non-extinction 
is equivalent to the initial condition (|Z^(0)|, |Z^(0)|) = (1,0), which says 
that the process starts with one skeleton particle, hence never dies out. For 
d > 1, let us define the range of Z^ as R^, and the range of Z^ as R^. 
Similarly, define 


M^{t) :=inf |r > 0 : R^{t) C Bj.(0)} , 
:=inf {r > 0 : R^{t) C B^(0)} . 
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Now since M{t) = max /*(0) = 0, (/*)'(!) — 1 = m, and 

/3m in the speed is invariant under how one describes the branching (see the 
remark after Proposition [1]), Proposition [1] immediately gives the desired 
lower bound: Ve > 0, 

P ^^ > y^2/3m — —)• 1 as t —>• oo. 

In order to obtain the desired upper bound, we find the mean matrix for 
the two-type continuous time branching process |Z| = {\Z^\, \Z‘^\). Let 

hj>0 

hj>0 

be the p.g.f.’s that determine the distribution of the offspring produced by 
particles of type 1 and type 2, respectively, where p^^\i,j)^ fc = 1,2 is the 
probability that a particle of type k gives i offspring of type 1 and j offspring 
of type 2. It follows that 


00/00 


f*{s) = Y 

i=0 \j=0 


and 

oo / oo \ 

j=0 \i=0 / 

where p^‘^\i,j) = 0 for i = 1, 2,... since a doomed particle never produces 
a skeleton particle. In order to find the mean matrix, we first calculate 
the expected number of offspring of type 2 for a particle of type 1, denoted 
by Li[Al(^)(l)|A^(0) = (1,0)], where the discrete-time process {N{n))nm 
represents the number of particles in generation n of \Z\. Let L denote the 
random variable for the number of offspring for a particle in the original 
GWP, L* of which are skeleton particles. By 0, Proposition 4.10], the 
joint p.g.f. oi L — L* and L is 






CONDITIONAL SPEED OE BBM AND RANDOM OBSTACLES 


18 


It follows that 

E[s^-^*t^‘'\L* > 1 ] 


ET=o,k=iPiL-L*=j,L* = k)sH^ 
P{L* > 1) 

f{qs + qt) - PiP = T L* = 0)s^ 
1 -q 

figs + qt) - ET=o PjL = j\L* = 0)P{L* 

1 - q 

f{qs + qt) - f{qs) 
l-q 


0)sJ' 


(32) 


where we have used (f^l) in the fourth equality. Setting s = t in ([3^ yields 
the p.g.f. for the total number of offspring for a skeleton particle: 


E[s^\L* > 1] = 

Finally, differentiating (f33l) with respect to s, setting s = 1 and then sub¬ 
tracting n gives the desired expectation 

E[N^^\l)\N{0) = (1,0)] = (34) 

It follows that the mean matrix [^, Section 5.7] for this two-type continuous 
time branching process is given by: 


A{t) := 





/3tm _ g 

gl3tk 



where Aij{t) is the expected population of type j at time t given that the 
process starts with a single particle of type i. 

The strategy now is to compare the two-type BBM Z at hand with a 
modified two-type BBM, say Z, prove the result for Z, and finally conclude 
that the result holds for Z by comparison. Let Z = (Z*^,2^) be a two-type 
BBM, P be the corresponding probability, and be the p.g.f.’s for 

the offspring of particles of type I and type 2 respectively, such that the 
p.g.f.’s for every 0 < s < 1 and every 0 < t < 1 satisfy 

=f^^\s,t) (35) 

OO / CO \ 

j=0 \i=0 J 


( 36 ) 
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which says that a particle of type 1 of Z and that of Z behave exactly alike, 
whereas a particle of type 2 of Z produces no offspring of type 1 similar to 
Z, but produces offspring of her own type just as a particle of type 1 of Z 
produces offspring of type 1, with mean > 1. The mean matrix for this 
modified process is 


Ait) := 




^/Stm 


Define similarly as before the range of Z^ as , and the range of Z^ as li. 
Further define 

M^{t) :=inf |r > 0 :^{t) C 5^(0)} , 

Ffit) :=inf |r > 0 :^(t) C 5^(0)} . 

By ([35]) and ([Mil , and since (/*)'(!) — 1 = m and /*(0) = 0, Proposition [T] 
implies that for every e > 0 


+ A = P + ^0 as t 


oo. 


Furthermore, since particles of type 1 of Z behave exactly the same way 
as those of Z, and particles of type 2 of Z produce offspring according 
to f* while particles of type 2 of Z produce offspring according to /, by 
Lemma [5] and Lemma [6l to complete the proof it suffices to show that 

P > \/2/?m + —)• 0 as t —)• oo for any e > 0. 

We proceed as in the proof of Proposition [H Let B = ^(^ 2/?m+£ )t(0)i 
5 > 0 satisfy (l27)l as before. Define'^^ = inf > 0 : |supp(Z^ (t)) n > l|, 
and estimate as before, 

p(^^ > + < P{\^it)\ > 

+ Pi4’B<t\ |^(t)| < 

By the fact that E[\^\ \ |Z(0)| = (1,0)] = |(m — k)j3te^^'^ and the Markov 
inequality, the first term on the right-hand side tends to zero. The second 
term on the right-hand side tends to zero by an extension of LemmaO Note 
that since each particle of type 2 of Z produces at least one particle of its 
own type (i.e., /(^)(s,0) = 0), we may extend Lemma [3] to apply to particles 
of type 2 of Z. This completes the proof of Theorem [TJ 
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5. Proof of Theorem [2] 

Part 3 of the theorem follows since r < oo almost surely for a non- 
supercritical process. For a detailed proof of part 3, see [l^, Theorem 1], 

Since A(0) = 0 gives a = 1 and P{£'^) = 1, conditioning the process on 
is the same as not conditioning it when A(0) = 0. Hence, we may extend 
part 2 to cover for the case A(0) = 0, and prove the first two parts of the 
theorem together. Note that the key ingredients in the proof are results 
concerning the decomposition of supercritical GW-processes and the speed 
of BBM, namely Lemma 0] and Theorem [TJ 

5.1. Proof of the lower bound 

Let m > 0. Fix /3, d, f and r], c. The scenario below yields the de¬ 
sired lower bound for the annealed survival probability conditioned on non¬ 
extinction, i.e., (E X P) {T > 

1. Recall that the first particle must be of infinite line of descent once the 
process is conditioned on non-extinction. Suppress the branching so 
that there is precisely one skeleton particle in the system up to time 
r]t. This is equivalent to requiring that whenever a skeleton particle 
branches within the period [0, ijt], it gives precisely one offspring of its 
own kind. Therefore, by the fourth part of Lemma [H this event has 
probability 

exp[—f3ar]t]. 

2. For d > 2, empty the two-sided cylinder (“tube”) Tt, defined by 

Ti = |x = {xi ,... ,Xd) e : |xi| < kt + h{t), < r{t) + /i(t)| , 

where k > c, and t i—)• r(t) and t i—)• h{t) are non-negative mappings 
that are picked such that lim^^oo= 0) = oo and 

h{t)/t = 0, limj^oo h(f) = oo. By the upper bounds above 
on r{t) and h{t), Lemma [T] gives that the probability to empty T^ is 
exp[o(f)] ford > 2. Ford = 1, empty the line T/ := {xi G M : |xi| < 2/i(t)}, 
which again costs a probability of exp[o(t)]. 

3. For d > 2, move the single skeleton particle (see the remark after 
Lemma ID during the period [0, rjt] so that it is at a specific site at 
distance ct + o{t) from the origin at time rjt (assume without loss of 
generality that the specific site has first coordinate ct + o{t) and all 
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other coordinates zero), and confine it to the smaller tube Tt, defined 
by 

ft = <x = {xi, ...,Xd) G : |xi| < kt, \Jxl + ... + x\ < r{t) I , 


up to time rjt. 

Let Pq be the d-dimensional Wiener measure, be the projection 
of the d-dimensional Brownian motion onto the first coordinate axis, 
and decompose the Brownian motion into an independent sum W = 
, where is hence implicitly defined. Define the events 


At = {ct < \W^t\ <ct + o{t)} , 

Bt = {\W}\ </ctV0<s<r/t}, 

Ct = < r{t) V 0 < s < rjt^ . 


Note that the event n P* n Ct is exactly the desired scenario for 

this part. The formula for the transition density for d-dimensional 

2 1—1 

Brownian motion gives Po{At) = exp[—^t + o{t)], and Lemma[7] gives 
Po(Pf) = exp[—-|- o{t)], so it follows that 


exp -—t + o{t) =Po{At) > Po{At n Bt) > Po{At) - Po(Pf 
2 rj 

= exp 


— exp 


= exp 

2 




which yields 


Po{At nBt) = exp 


--t + oW 


(37) 


Since r{t) —)■ oo, we have Po(C't) = exp[o(t)]. Combining this with 
(|37|) and using the independence of and we find 


Po(^t n Pt n Ct) = exp 


--« + o(() 


For d = 1, note that the function in ([5|) is minimized when c = 0 
(see the proof of Theorem [3|) . Confine the single skeleton particle to 
the one-dimensional tube T/ := {xi G M : [xil < h{t)} up to time rjt. 
Since the length of the tube tends to infinity as t tends to infinity, the 
probability of this event is exp[o(t)]. 
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4. Confine all the doomed particles that are created within the period 
[OjT/t] to the larger tubes Tt and T/ for d > 2 and d = 1, respectively, 
up to time t. Since the number of occurrences of branching up to time 
t along a single ancestral line is a Poisson process with mean /It, and a 
skeleton particle on average produces {fi — K)q/{l — q) doomed particles 
every time it branches (see the proof of Theorem [T|) , it follows that at 
most [f3rjt{fj.—K)q/ (1—g)J doomed particles are produced along the sin¬ 
gle skeletal line up to time rjt with probability exp[o(t)]. Let the radi0 
of the subtrees initiated by these — K)q/{l — q)\ =: n{t) doomed 

particles be pi,p 2 ,-- ■,Pn{t)- If K{t) := max{pi,p 2 , • • -^Pnit)}, then 

PiK{t) < h{t)) = P{pi < (38) 


Now, since a doomed subtree is almost surely finite, lim^^oo h{t) = 
oo implies that P{pi < h{t)) 1, which in turn implies that the 

probability in (l38|) is exp[o(f)]. We have already confined the single 
skeleton particle to the smaller tube Tt for d > 2 (P/ for d = 1). Since 
each doomed particle that is produced along the single skeletal line 
within the period [0, rjt] must be produced within the smaller tube 
Tt{T^), and the dimensions of the larger tube Tt{T^) all exceed the 
corresponding dimensions of the smaller tube by h{t), we conclude 
that all of the doomed particles that are created within the period 
[ 0 ,r 7 t] stay inside the larger tube with probability exp[o(t)]. 

5. Empty a y/2flm{l — r])t-haH around the position of the single skeleton 
particle at time pt. By Lemma [H this event has probability 


exp 


-^d<i(\/2/3m(l - p),c)t + o{t) 


6. Require the branching system initiated by the single skeleton particle 
present at time pt to stay inside the y/2^rn{l — p)t-hall during the re¬ 
maining time {l — p)t. When m > 0, by Proposition [T] and Theorem [H 
the probability of this event conditioned on non-extinction is exp[o(t)]. 

Since the motion, branching and trap formation mechanisms are indepen¬ 
dent of each other, minimizing the cost of all these events over the parame¬ 
ters p and c provides us with the desired lower estimate for (IExP)(r > t\£^). 
Note that the survival scenario described above is composed of three large 


^By radius we mean the radius of the smallest ball centered at the root of the subtree, 
containing the (finite) range of the subtree. 
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deviation events: suppressing the branching for a linear time, moving the 
Brownian particle to a linear distance, and emptying a ball with linear ra¬ 
dius. The remaining three components, namely parts two, four and six are 
not large deviation events and do not give exponential costs. 


5.2. Proof of the upper bound 

Fix (3, d, f. Let 0 < e < 1. Recall that \Z(t)\ is the number of particles 
in the system at time t, |Z*(t)| of which are skeleton particles. For t > 1, 
define 

= sup |r? G [0,1] : \Z*{r]t)\ < | . (39) 

Before we proceed, we prove the following lemma: 


Lemma 8. 



= exp[—/3a—t + o{t)]. 


Proof. Observe that < Vt} = {\Z*{^t)\ < By part 4 of LemmalU 

we know that the process \Z*\ conditioned on non-extinction of |Z| is equal 
in law to a continuous time branching process with rate /3(1 — f'{q)) = /3a 
and offspring distribution A, where A(0) = A(l) = 0. On the event 
since the first particle is of infinite line of descent and t > 1, the inequality 
< 'nt\£^) > exp[—,5a^t -|- o{t)] follows trivially. 

Now consider a strictly dyadic branching Brownian motion Z. We have 
P(\ Z(t) \ = k) = e~^^{l — for every k G N and for every t > 0, 

see [l^. It follows that P{\Z{t)\ > k) = {1 — Using the binomial 

expansion (1 — x)” = 1 — nx -|- — ... (—l)"'x”, we obtain 

p(|z(t)i < 

=e-^* + • • • + 

^ exp[-/3t -b o{t)]. 


The result follows by comparison with the strictly dyadic case, since A(0) = 

A(1)=0. □ 


Henceforth, all probabilities written should be understood as conditioned 
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on Now, for every n S N, 


(E X P){T > t) 


n—1 


2=0 

n—1 

< ^exp 
i=0 


i + 1 


xP)({T>t}n^-<r,t< + (E X P) ({T >t}n ivt = 1}) 

(E X {T > f} + exp[-/3t + o(f)], (40) 


—fla—t + oit) 
n 


where we use Lemma [H and introduce the conditional probabilities 




■)=P 


i i + 1 

-<'nt< - , ^ = 0,1,... ,n - 1. 

n n 


Recall that {rft < (i + l)/n} = + l)/n)| > Let A^^’'^\ i = 

0,1,... , n — 1 to be the event that among the + l)/n)| skeleton 

particles alive at time t{i + l)/n-, there are < particles such that the 

ball with radius 

^Rn) (1 _ ^)^2/3m(l - (41) 

around the particle is non-empty (i.e., contains a point from oj). Estimate 

(E X P^(*’"))(T > t) < (E X p/*’”^)(^f’”^) 

+(E X p^(*’”))(r > t I [A^t'^^^f)- (42) 


On the event there are > balls containing a trap, and by 

Proposition [1] and (1411) , the BBM emanating from the center of each ball 
exists this ball in the remaining time (1 — (i + l)/n)t with a probability 
tending to 1 as t —^ oo. By radial symmetry, the BBM hits a trap inside 
this ball with a probability when exiting, where Ci > 0 is a 

constant depending on the geometry. Hence, the second term on the right- 
hand side of (l42]) is bounded above by 


1 - 


Cl 


1 


Ip!”"’]"-' 


< exp[—C' 2 t 


l + dl 


uniformly in all parameters, where (72 > 0 is another constant. Hence, the 
second term on the right-hand side of (|42l) is superexponentially small (SES). 
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Now consider the first term on the right-hand side of ()42p . Since gdi^t, bt) = 
tgd{r, b) for r,b > 0, the cost of clearing a ball with radius linear in t is ex¬ 
ponentially small in t as t —)• oo. Then, for large t, we have 

(E X I \Z*{t{i + l)/n)| = +j) 

> (E X I + l)/n)| = + (j + 1)) 

for every j G {1, 2,...}. It follows that 

(E X < (E X I \Z*{t{i + l)/n)\ = -b 1). 

(43) 

We continue to estimate from above as follows. Each skeleton particle alive 
at time is at a random point, whose spatial distribution is identical 

to that of W{t{i + I)/n), where W denotes the standard Brownian motion. 
Let xq be the center of the empty ball at time t{i -|- l)/n that is closest to 
the origin. Let 0 < c < oo, 5 > 0, and define Bj, j = 1,..., -|- 1 to 

be the event that jth skeleton particle at time t{i -|- l)/n is at a distance 
ct < r < {c + 5)t from the origin and that the pj^’^^-ball around it is empty. 
Then, 

n {ct < |xo| < (c + <5)t} I \Z*{t{i + l)/n)\ = + 1) 

<(ExP/'’”))(BiU...UBL,.+ej+i) 

+ 1)(E X 


=(Lt'^+"J +l)exp 
X exp 


-t + o{t) 


2 {i + l)/n 

Igd f(l -e)\/2^(m - 1)(1 




),c ) t + 0{6)t + o{t) 


(44) 


where we have used the independence of BBM and trap mechanisms, and 
Lemma [T] in passing to the last equality. Indeed, on the event condi¬ 
tioned on \Z*{t{i + 1)/n)\ = -|- I, there is at least 1 skeleton particle 
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with an empty ball around it. From (j42p . (l43D and (I44p . we obtain 


(]ExP/*’''^)(r > t) 


n—1 


< ^(E X pf I ^ < |xo| < ^ [ 

j=0 

+ (E X 

n—1 

, l‘M~ ! Tl~ 

t + o{t) 




n A 


{i,n) 


< 


+l)exp 


({|xo| > v^t} n + SES 


j=o 


{i + l)/n 


-Igd f (1 - e)y/2Pm{l - i + 0{l/n)t + o{t) 


X exp 


+ exp[—+ o{t)] + SES. 


n n 


(45) 


Here, the SES comes from the second term on the right-hand side of (I42p . 
Also, in passing to the last inequality, we have used that the probability of 
the event {|xo| > \/2^t} is bounded above by the probability that a single 
Brownian particle is at a distance > ^/^t at time t, which is exp[—/3t-\-o{t)]. 

Substituting (HSl) into (HOl) . and optimizing over i,j G {0,1,... , n — 1} 
gives 


limsup-log(E X P)(r > t) 

t—>-CXD t 


< — min 


Pai Pf/n^ 

n (i + l)/n 


+ ^9(1 



— e)^j2l3m{l 




Now let rj = i/n, c = j/ny/^, let n —oo, use the continuity of the 
functional form from which the minimum is taken, and finally let e —)■ 0 to 
obtain the desired upper bound, 

limsup- log(E X P){T > t) < —I{1, /, /3, d). 

t—>-oo t 


Note that c > y/^ cannot minimize 

{/3af? + ^ + l9d{\/Wm{l - rj), c)| (46) 

since this makes (I46p greater than (3, which can be improved by setting c = 0 
and T] = 1. 
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6. Proof of Theorem [3] 


Let 

c 2 ,_ 

Gdiv, c) = liar} + ^ + 2/3m(l - t/), c) (47) 

so that 

min Gd{v,c). (48) 

r;e[0,l],ce[0,V2?] 

To see the existence of the minimizers r]*, c* of (|48p . note that Gd is lower 
semicontinuous on the compact subset [0,1] x [0, \/2p\ of since 


lim 

(r,,c)^(0,0) 


I3ar] + Igd (^V2/3m(l - t]), = lsd\/2^m 


and g > 0 . 

We refer the reader to (7| for the proof of their uniqueness when I ^ Icr- 
Consider d = 1. Since gi{r, b) = 2r, the minimum over c in (I48p is taken 
at c* = 0 , so that (|48]l reduces to 

/(/,/,/?,d) = min |flag + 2/y^2/3m(l — r?)| 

»76[0,1] t J 

= min Ir/(/3a — 2l^2j3m) + 2l\j2fim \ . 

?7e[o,i] I J 

This proves ® and ([T|), and identihes Icr as the solution to /3a = 2ly/2/3m. 
Now consider d> 2. We have 

Gdiv, c) - Gd{0, 0 ) = Par] + — 

2g 

+ 45 d(\/ 2 / 3 m(l - g),c) - 5 rf(v^ 2 / 3 m, 0)], (49) 

where the last term on the right-hand side is less than or equal to zero, 
with equality if and only if (77, c) = ( 0 , 0 ) (this follows from the definiton 
of gd)- Suppose that {g,c) = (0,0) is a minimizer when I = Iq. Then, by 
uniqueness of minimizers, the right-hand side of (1491) is strictly positive for 
all {g,c) 7 ^ ( 0 , 0 ) when I = Iq. As the last term in the right-hand side of 
(f4^ is strictly negative for all {g,c) / ( 0 , 0 ), it follows that for all I < Iq, 
the right-hand side of (14^ is zero when (r/, c) = ( 0 , 0 ) and strictly positive 
otherwise. Therefore, there must exist Icr S (0, 00 ) such that 

1 . (? 7 , c) = ( 0 , 0 ) is the minimizer when I < Icr, 

2. (? 7 , c) = ( 0 , 0 ) is not the minimizer when I > Icr- 
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Henceforth, we will take this as the definition of Icr for d > 2. This proves 
the first parts of ([8|) and ([9|). Note that Icr S (OjCo), because the last term 
on the right-hand side of decreases without bound as ^ oo and tends 
to zero as Z —)• 0 for fixed {rj, c) / (0,0). 

Now let d > 2 and I > lev We know that {ri,c) = (0,0) is not a 
minimizer. The combination rj* = 0, c* > 0 is not possible because G(0, c) = 
oo for all c > 0. The combination r]* > 0, c* = 0 is ruled out as well because 
G{r],0) takes its minimum either at ?? = 0 or r/ = 1, so ry* >0 would imply 
that rj* = 1. However, rj* = 1 can be excluded via [3, Lemma 4.1], Also, 
c* = is not possible since this yields > 13, whereas {t],c) = (1,0) is 
not a minimizer yet yields G^ < (3 since a £ (0,1]. Hence, we conclude that 
the minimizers for d> 2 and I > Icr appear in the interior of [0,1] x [0, y/2(3\. 
This proves the second parts of ([8|) and dH). 

In the rest of the proof, we find Icr for d> 2. For i? > 0, let 


gd{R) = [ 

JE 


dx 


Br( 0 ) \x + e 


d-i ■ 


(50) 


Then, we may write (1471) as 


Gdid, c) = Par] + — + IcQd 
If] 


V2/3m(l - Tj) 


(51) 


The stationary points are the solutions of the equations = 0 and = 
0, which yield 

0 =/3a — /3mv‘^ — l\/2f3mg’^{u) 

0 =y/2l3mv + l[gd(,u) - ug'diu)] (52) 


upon setting 

y/2(3m(l — n) c 

u =-, V = ■ 

c y/2pmg 

Eliminating g'^ in (|52p gives 


1 ^ 9 d(u) = -v + -^u{a - mv'^). 
y/2f3m 2m 


(53) 


(54) 


It follows by (fKTI) . (|^ and (|Ml) that for d > 2 and and I > Icr, 


I{l,f,f3,d)=Gd{u*,v*) = l3{a-mv*^). (55) 


Since the minimizer [u* ,v*) must satisfy ([54]) (note that there may be more 
than one pair {u,v) that satishes (l5H) L we conclude by (1551) that v* is the 
maximal value of u > 0 on the curve in the {u, u)-plane given by (|54p . 
Recall the following lemma from [^, where a proof is given as well. 
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Lemma 9. i? i—)• g'^iR) is strictly increasing on (0,1), infinite at 1, and 
strictly decreasing on (l,oo). 

Using ([3|), we may write (IMl) as 


gd{u)l 

2sdml*r 


—V -\ - u{a — mv^). 

2m 


(56) 


Since gd{u) ~ SdU as u —>■ oo by (|50l) . we obtain from Lemma 0 that the 
left-hand side of (fSGl) is strictly convex on (0,1), has infinite slope at 1, is 
strictly concave on (l,oo), and has limiting derivative l/{2ml%) (see Figure 
1 ). 



Figure 1. Qualitative plot of: (1) u !->■ —v + -^u{a — mv"^). 

The dotted line is u 

2m i* 


Firstly, see that Icr < al*^. Indeed, assume the contrary. Then, there 
exists I G {allr, fir), which implies that {g, c) = (0,0) is the unique minimizer 
for Gd when I = 1. However, Gd(0,0) = fil/fir > fioi = Grf(l,0), which 
contradicts {g,c) = (0,0) being the minimizer. This implies that fir < afir- 
We claim that fir can be identified by the following formula 


^0 


sup 


I > 0 


I 1 

fir* 2sdm 


gd{u) > 



1 I 

2m fir* 


Vu € (0, oo) 


(57) 

In order to prove our claim, i.e., fi = fir, we need to show that (r/, c) = 
(0,0) is the minimizer for Gd when I < fi, and (??, c) = (0,0) is not the 
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minimizer when I > Iq. Equivalently, as G(i(0,0) = fil/l*cr-, we need to show 
via (f55]l that for I > Iq there exists a v > ^Jljmia — l/lcr*) such that v 
satisfies (1561) for some u > 0, and that for I < Iq there is no such v. As we 
know that Icr > 0 and (r/, c) = (0,0) is the unique minimizer when I < lc.r^ 
there must exist an /i >0 such that the inequality in (j57l) holds for all 
u € (0, oo) when I = h, hence 0 < < Iq. Now take any 0 < I 2 < h, and 
see that the inequality in (l57)) holds for all u G (0,oo) also for I = 12- This 
implies that if I < lo, then the curve on the left-hand side of (|56p and the line 
on the right-hand side do not touch when v = ^Jljraia — l/lcr*)- As both 
terms on the right-hand side of (1561) decrease as v increases and the left-hand 
side of (IM]l does not depend on v, the curve on the left-hand side of (IMI) 
and the line on the right-hand side do not touch for v > \J\/m{a — l/lcr*) 
as well when I < Iq. 

Now consider I > Iq. By definition of Iq, there exists u G (0,oo), say u, 
such that 


h 1 ( \__ \}_( ^0 ^ 1 Iq 

Icr* y m \ Icr* ) 2m Icr* 


(58) 


Note that since gd{u) < s^u for all u G (0, 00 ), (l58l) implies that we have the 
strict inequality Iq < Now take any li G (Zo,a^cr)i observe that at 
u = u, the right-hand side of (1581) becomes greater than the left-hand side 
if we replace Iq by li. Also, since gd{u) ~ SdU, the left-hand side of (1581) 
must be greater than the right-hand side for large u if we replace Iq by li. 
As the functions on both sides of ([581) are continuous with respect to u, it 
follows by the intermediate value theorem that there exists u G (0,oo), say 
u, such that (1581) holds, with Iq replaced by li. This in turn implies that 
Gd{u,v) = Gd{u, Y^l/m(a — h/lcr*)) = l^h/llr^ since the minimizer is 
unique, this shows that {rj, c) = (0, 0) is not the minimizer for Iq < li < 

As we have shown that ( 77 , c) = (0,0) is not the minimizer when I G {Iq, Icr), 
we conclude by definition of Icr that ( 77 , c) = (0,0) is not the minimizer for 
Gd for all I > Iq. Hence, Iq = Icr- 
Now put 

gd{u) = SdU - gd{u) Aid = 77^ max gd- 

ISf^TTl liE(oD) 


Then, (|57)) reads 



IN 


Icr = sup < / > 0 : 


(59) 
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This implies that 



which completes the proof of (fT0|) and (fTT]l . 

We finally show that u* G (0, 1) when I > Icr- Note that for a fixed 
V, if the line on the right-hand side of (1561) cuts the curve on the left-hand 
side at more than one point, then by (I55p . this v cannot be v* since the 
minimizer is unique. Hence, we are looking for a v value such that (1561) is 
satisfied for exactly one u value, and this pair is {u*,v*). This implies that 
when V = V* , the line on the right-hand side of (1561) is tangent to the curve 
on the left-hand side at u = u*. Since v* > ^tjmia — l/lcr*), the slope 
of the line on the right-hand side of (1561) is < l/{2m){l/l*r). Then, since 
the left-hand side of ([56]) has infinite slope at u = 1, is concave on (l,oo), 
and decreases asymptotically to l/{2m){l/l*^), it follows that u* G (0,1). 
By (l53]) . we conclude that c* > y/2^m{l — Vj*). This completes the proof of 
Theorem [3l 
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